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ABSTRACT 

We explore the halo-to-halo variation of dark matter substructure in galaxy-sized dark matter 
hales, focusing on its implications for strongly gravitational lensed systems. We find that the median 
value for projected substructure mass fractions within projected radii of 3% of the host halo virial 
radius is approximately fsnh ~ 0.25%, but that the variance is large with a 95-percentile range of 
< /sub < 1%. We quantify possible effects of substructure on quadruply- imaged lens systems using 
the cusp relation and the simple statistic, -Rcusp- We estimate that the probability of obtaining the 
large values of the i?cusp which have been observed from substructure effects is roughly ^ 10"'^ to 
^ 10^^. We consider a variety of possible correlations between host halo properties and substructure 
properties in order to probe possible sample biases. In particular, low-concentration host dark matter 
halos have more large substructures and give rise to large values of -Rcusp more often. However, there 
is no known observational bias that would drive observed quadruply-imaged quasars to be produced 
by low-concentration lens halos. Finally, we show that the substructure mass fraction is a relatively 
reliable predictor of the value of -Rcusp- 

Subject headings: galaxies: halos - gravitational lensing: strong - theory: dark matter 



1. INTRODUCTION 

In the standard Cold Dark Matter (ACDM) cosmology, 
structure forms hierarchically. Small dark matter halos 
form first and merge to form larger objects. The process 
of merging and dynamical evolution leaves remnants of 
smaller halos (called subhalos) within larger host halos. 
We refer to the population of subhalos as substructure. 
The presence of substructure in galaxy-sized dark matter 
halos is an important test of cold dark matter. While the 
distribution and number of galaxies and clusters of galax- 
ies is in agreement with the generic predictions of CDM, 
at sub-galactic scales numerical simulations predict an 
abundance of dark matter substructure which has not 
been detected by the presence of baryonic matter (stars 
& gas) , a situation k nown as the "missing satellite s prob 



et al. 



lem" (Klypin et al. 1999; Moore et al. 



1999: Diemand 

i ■ 



jUU7| |apringel et al. 2008'). Detectmg dark mat 



ter subhalos and putting constraints on the substructure 
mass function is a key step towa rd the identification o f 
the dark matter particle (see, e.g., Moustakas et al. 2009 ). 

One of the most promising avenues to measuring sub- 
structure in galaxies is via strong gravitational lensing. 
As light from a distant object passes an intervening 
galaxy, it may be bent sufficiently for multiple images 
of the source to be observed. Strong lens systems with 
quasar sources and galaxy-sized lenses have been shown 
to be sensitive to satellite galaxv-sized dark matter sub- 
halos in the lensing hal(f]|Mao fc Schneider|l998||Metcalf 
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2002 |Chiba|[2002| [Kochanek fc Dalal||2004| ). Substruc- 
ture m the lensmg halo may be observed as perturbations 
to the image magn ifications (flux ratio anomalies) , im- 



age p ositions (e.g., Metcalf & Madau 2001 Chen et al. 
2007 1, a nd/or time delays in system s with time variable 
sources (Keeton & Moustakas| 2009|). 



The most well-explored approach to substructure lens- 
ing is the study of anomalous flux ratios in multiply- 
imaged systems. Anomalous flux ratios are sensitive to 
the projected mass fraction in substructure, /sub, at pro- 
jected radii where lensed images are present, ~3% of the 
projected virial radius of the lens halo. Results , how- 
have not been definitive. IDalal & Kochanek (2002 1 



ever. 



find that /sub = 2% with 0.6% < /sub < 7% withih"Wo 
confidence intervals. This is margina lly consistent with 
the re sults of numerical simulations. IMetcalf fc Amaral 
( 2010 1 also find a result which is generally consistent with 



numerical simulations. 

Several studies have shown that /s ub is small (< 1%) 
but wi th large variance. For example, |Zentner fc Bullock] 
(20031 used semi-analytic simulations of a large num- 
ber or halos to argue that substructure mass fractions 
at projected radii of 3% of halo virial radii should be 
■fsub ~ 0.6% on aver age and /sub ^ 1% for 95% of ha- 

< 0.5 % using 1 2 halos 



los. Mao et al. 



(|2004) found /sub 
id cluster masses. 



at galaxy, g roup ani 

found .01% < /sub < -7%, using 6 halos f rom the Aquar 



Xu et al. (20091 



ius simulation. Using 26 simulated halos, Zentner (20061 
showed that /sub < 1% along most lines of sight, but that 
there is a large variance due, in part, to the orientation 
of the halo with respect to the line of sight. 

Other studies have suggested that lensing observa- 
tions are inconsistent with expectations from numeri- 
cal simulations. Using simple models for substructure. 



Chen (20091 found that observations support the pres- 



more substructure than expected from simula- 



components within that dark matter halo. 
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tions. Maccio et al. (2006) and Amara et al. (|2006l 
each compared observations to a single simulate 



halo 



and found more and greater magnification perturbations 
than can be accounted for in their simulated hal os. In 
addition, the results of Dalai & Kochanck (2002) and 



Metcalf & Amara ( 2010 ) (which are consistent with sim- 
ulations) use substructure models that make it possible 
that their results underestimate the value of /sub- 

Future, more definitive, results will require better char- 
acterizations of sources of error, more observational data, 
and better predictions for CDM substructure. In this 
work we focus on this last element. More specifically, as 
the initial density field is unknown, predictions for sub- 
structure populations in numerical simulations are spe- 
cific to a particular realization of the simulated halo. 
Halo-to-halo variation, which should be present, is es- 
sential in determining the expected variance about each 
measurement. Quantifying this variance can allow more 
accurate error assessments and potentially help bring ob- 
servational results into better agreement with theoretical 
predictions. In addition, understanding the variance of 
the substructure population may provide valuable insight 
into previously unexplored biases in the lensing samples. 

In this work, we characterize the halo-to-halo variation 
of substructure populations in galaxy-sized dark matter 
halos and investigate how substructure properties cor- 
relate with other properties of dark matter halos. We 
study the effects of halo-to-halo variation to mock lensing 
observations, focusing on a simple statistic that quanti- 
fies violations of the cusp relation. We employ a semi- 
analytic method to generate a large number of substruc- 
ture populations. We describe this method in Sec.j2]and 
describe basic properties of halo substructure in Sec. [Sj 
We use our substructure models to generate strong lens- 
ing simulations in Sec.|4] We analyze our simulations and 
compare the results to observations in Sec. [5j we discuss 
caveats in Sec.|6j and we conclude in Sec. [7] Throughout 
this paper, we work in a flat cosmological model with 
Q,dm = 0.228, fife = 0.0227, h = 0.71, and erg = 0.81 as 
favored by the five-year Wilkinson M icrowave Anisotropy 
Probe results ( Komatsu et al.|[2009 l. 



2. SEMI- ANALYTIC MODEL 



In order to generate a large statistical sample of host 



halos , we use a semi-analytic technique |Zentner et al. 
2005 1 . This technique allows us to approximate the hi 



erarchical assembly of a dark matter halo analytically, 
and thus avoid the prohibitive computational cost asso- 
ciated with nu merous, high resolution num erical simula- 
tions (see also [Taylor fc Babul|[200l| |2004[ ). 

The semi-analytic approach employs a variation of the 
extended Press-Schechter formalism to generate in erger 
histories of host dark matter halos (Bond et al. 1991; 
Lacey fc Coie| [T9931 [Somerville fc Kolatt||l999, ./ent-. 



ner 



2007[ ) aiig an approximate analytic approach to the 



dynamical evolution of subhalos in an evolvi ng host 



halo. The details of this method are given in Zentner 
et al. ( 2005 1 , including a number of comparisons to self- 



consistent numerical simulations. This technique has 
been tested for subhalos of mass M > 10"'' Mh, where 
Mh is the mass of the host halo, and produces subhalo 
abundances and spatial distributions in good agreement 
with A'^-body simulations. Though this approach is ap- 
proximate, it allows us to compute a large number of 



realizations of a host halo at a subhalo resolution which 
is inaccessible at present with state-of-the-art numeri- 
cal simulations. In addition, this procedure can be used 
to produce physically moti vated extrapolation s t o nu 
meric al simulations (see Koushiappas & Zentner" '2006 
,Koushiappas et al. 2004^ ,Koushiappas et al. _2010). 



We generate 200 realizations ot a host halo of mass 
Mh = 1.26 X lO^2/i"^M0. We track the accretion and 
dynamical evolution of subhalos of mass greater than 
10~^Mh. The mean properties of substructure in the 
sample of realizations is in agreement with numerical 
simulations (see Fig. 1 in [Koushiappas et al.||2010[ ). The 
average number of subhalos with Vmax > 4 km/s within 
the inner 200 kpc of the halo is 2481 with a 68 percentile 
range of [1964 — 3007], which is consistent with the 2469 
subhalos found within the same radius of a numerical 
simulation of a similar-size halo 



(Kuhlen et al. 2008 



In addition, the scatter in the number ot subhalos from 
halo-to-halo variation is consistent with a Poisson scatter 
added in quadrature to an intrinsic scatter of ab out 20% 
(as derived from re cent numerical simulations; Boylan- 
Kolchin et al][2010[ ) 



Flach realization represents a possible subhalo popu- 
lation within the host halo. The distribution, structure 
and properties of the subhalo population differ between 
realizations due to the statistical nature of the initial 
density field and the complexity of any individual sub- 
halo's orbit. The output of each realization contains the 
concentration of the host and a list of all subhalos, in- 
cluding their evolved structural parameters - total bound 
mass, scale radius, a.nd tidal radius - an d their locations 
within the host (see Zentner et al.|2005 ). The individual 
merger history of each host influences the abundance and 
properties of its substructure population. 

3. SUBSTRUCTURE STATISTICS 

We construct a two-dimensional projection of the sub- 
halo distribution by projecting the three-dimensional 
density along an arbitrary direction. The three- 
dimensional density profile of each subha fo is described 
by a Navarro-Frenk- White profile (NFW; Navarro et al. 
1997 1, p{r) — ps/[x{l + x)^], where x = r/vs, rg is the 
scale radius, and r is the three-dimensional distance from 
the halo center. Each subhalo has a projected distance 
from the center of the halo, which we call R. For sim- 
plicity we define a dimensionless projected distance for 
each subhalo as R = R/R^, where i?h is the virial radius 
of the host halo. This two-dimensional projection of the 
subhalo population is used in the mock lensing simula- 
tions we describe below. We also utilize a dimensionless 
subhalo rnass M defined such that the mass of a subhalo 
is M = MMii, where Mh is the mass of the host halo. 

We can glean information about the variation in the 
distribution of substructure by calculating the subhalo 
mass fraction within a projected radius. The projected 
mass fraction within the inner R fraction of the virial 
radius of the host in a given subhalo mass bin is 

/sub(< R) = ~ / ~, E < ^) 

Mh[< R) J 



where 



Mhi< R) 



-^(c) Jo ^ JO 



f[l + f] 



di dR! (2) 
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Fig. 1. — The projected substructure mass fraction as a function 
of radius (R = R/R\i) for sublialos with mass greater than 10^^ 
times the mass of the host halo. The median in the sample of 
200 realizations is shown as a point, the 68 percentile is shown 
with narrow errorbars, while the 95 percentile is shown with wide 
errorbars. The large spread of /s^b at large radii is an outcome 
of the distribution of merger histories of the host halos. The top 
panel shows the projected substructure mass fraction at radii R = 
[0 — 1], while the bottom panel shows the projected substructure 
mass fraction in the inner regions of the halo, which is relevant to 
lensing studies. 

is the projected mass within the inner R fraction of the 
virial mass of the host halo. In Eq. [T] the sum over 
subhalo masses is over ah subhalos i whose projected 
radial position Ri is less than the projected host position 
i?, independent of the projected size of the subhalo. In 
Eq. [zj f = cr/i?ii = V R'"^ + the upper integral in the 

line of sight element is £max — — R'^ , and J-{c) — 
ln(l + c) — c/(l + c), where c is the concentration of the 
host halo. Note that values of /sub > 1 are possible. This 
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Fig. 2. — Top: The projected substructure mass fraction as a 
function of radius {R = R/R^^) for subhalos with mass greater 
than 10^^ times the mass of the host halo. Host halos with con- 
centrations in the 68 percentile of the mean have substructure mass 
fractions which are of order ~ 10%, while lower and higher con- 
centration hosts have a substructure mass fraction which is higher 
and lower, respectively. See text for details. Bottom: The mass 
fraction within the projected inner 3% radius of objects in a given 
mass bin for hosts with low (left panel), mean (middle panel) and 
high (right panel) concentration. In both figures, the y-axis error- 
bars represent the same percentile range as in Fig.[l] Points depict 
the median of the distribution. Note that as the concentration of 
the host increases, the fraction of mass contributing to the inner 
3% of the projected radius becomes negligible. As a result, only 
hosts with low concentration are likely to have rich substructure in 
the projected inner 3% of their radius. 

may be the result of two effects: Mh(< R) from Eq. [l] 
only refers to the smooth mass of the host halo and not 
the total mass, and we use the full subhalo mass even 
when most of the subhalo's bound mass falls outside of 
the relevant region. 

In the top panel of Fig. [T] we show the projected cu- 
mulative substructure mass fraction /sub, in the complete 
sample of 200 realizations as a function of radius. The 
median value of /sub at i? = 1 is 9% for subhalos with 
mass greater than 10~^Mh, with a 95 percentile in the 
range [2 — 17]%. Note that the spread about the median 
increases as R increases. This is due to the fact that the 
outer regions of halos are dominated by recently accreted 
subhalos, while the subhalo population of the inner re- 
gions is old and highly evolved. The recent accretion of 
subhalos is an outcome of the individual accretion his- 
tory of each host halo, and the large spread about the 
median in the outer regions mirrors the range of recent 
accretion histories. 

Relevant to lensing studies is the projected substruc- 
ture mass fraction within the inner 3% of the virial radius 
of the host halo. The bottom panel in Fig. [T] shows /sub 
in the inner regions. At R = 0.03, the median value of 
/sub is 0.25% with a 95 percentile in the range [0 — 1]%, 
meaning that only 2.5% of the realizations have a pro- 
jected substructure mass fraction greater than 1% in the 
inner 3% of the virial mass of the host halos. 

In order to explore the relationship between the host 
halo properties and the substructure lensing effects, we 
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distribute realizations into three groups. In the first 
group, which we caU Co^we include host halos with a 
concentration c (c.f., Eq. [2]) within the 68% range of the 
mean concentration of the complete sample, (c) = 9.7. 
This includes host halos with concentrations in the range 
6.7 < c < 13.4. We then collect the upper 16% of the con- 
centration range (c > 13.4) in group C+, and the lower 
16% range (c < 6.7) in group C_. As they form earlier, 
subhalos present in high concentration hosts are expected 
to be fewer and with more concentrated subhalo density 
profiles relative to subhalos present in low concentration 
hosts. The inverse relation between substructure abun- 
dance and halo concentration has also been seen in the 
simulations of ? and ?, where they measure the abun- 
dance of subhalos around dark matter halos of ~ lO^^M© 
and larger. Grouping the hosts (and their substructure 
populations) allows us to investigate the effects of host 
concentration on the expected lensing signal from its sub- 
structure. 

The cumulative substructure mass fraction for halos 
with mass greater than 10~^ the mass of the host is also 
dependent on the concentration of the host halo. In the 
top panel of Fig. [2] we show the median and 68 & 95 
percentiles of the distribution of mass fractions found in 
our sample. Host halos with concentration in the 68 per- 
centile about the median concentration of all realizations 
(middle figure) have a substructure mass fraction which 
is of order ~ 10%, and a 68 percentile of ~ [5 — 20]% 
(the 95 percentile range is ~ [3 — 35]%). Higher concen- 
tration hosts have a smaller fraction (of order ~ 6% and 
a 68 percentile range of ^ [3 — 10]%, and a 95 percentile 
range of ~ [2 — 25]%), while low concentration hosts have 
a higher substructure mass fraction (median of ^ 22%, 
and a 68 percentile range of [9 — 50] % and a 95 percentile 
of~ [4—55]%). Systems of low concentration (and subse- 
quently large values of /sub) are systems which typically 
contain recent large mergers. Such systems may appear 
as interacting galaxy pairs or groups in observations. 

In the bottom panel of Fig. [2] we show the projected 
mass fraction within the inner 3% of the projected virial 
radius (most relevant to lensing studies), for a number 
of mass binsH] We find that the abundance of subha- 
los with masses greater than ~ few times 10~^ the mass 
of the host is negligibly small: even though the upper 
95 percentile of the mass fraction increases with mass 
in the panels, the median and the upper 68 percentile 
drop to zero (as does the number abundance) as mass 
increases. Therefore, large subhalos contribute little to 
projected mass fractions on average, although, in the rare 
cases in which a large subhalo is present, the projected 
mass fraction is correspondingly large. This deficit of 
large subhalos is independent of the host halo concen- 
tration. However, on smaller scales (masses of order 

M ~ lO"'^ — 10~^), the effects of the host concentra- 
tion are prominent. For low-concentration host halos, the 
median mass fraction is /sub ~ 0.3% with an upper 68 
percentile of /sub ~ 0.8% (95 percentile is /sub ~ 1-3%), 
within the projected inner 3% of the virial radius. In 
contrast, the high-concentration host halos have a me- 
dian of /sub ~ 0.035%, and with an upper 95 percentile 
of less than 0.2%. 

''' This particular choice of mass bins is for illustrative purposes 
only. 



The origin of these differences between low- and high- 
concentration halos is the same as alluded to earlier. 
High-concentration halos generally form earlier than low- 
concentration counterparts, undergoing the majority of 
their mergers at higher redshift. As a consequence, 
the merging subhalos evolve for a relatively longer pe- 
riod within the dense environments of their hosts in 
high-concentration halos. This allows more time for 
orbital decay and mass loss. The increased propen- 
sity to lose mass in high-concentration halo is com- 
pounded by the fact that tidal forces are moderately 
stronger in high-concentration host halos at fixed, to- 
tal bound mass. These effects, working together, ren- 
der the survival rate of mass bound to halo substruc- 
ture within high-concentration hosts smaller than in low- 
concentration hosts. The opposite is true of substructure 
in low-concentration host halos, so low-concentration 
host halos tend to have higher substructure mass frac- 
tions. These concentration-dependent eff ects are consis- 
tent w ith those found and discu ssed in IZentner et al.l 
((20051) andlWechsler et all (l2006l). 



0.02 



V 0.01 





0.03 
iS 0.02 

V 

J 0.01 


0.06 

0.05 
S 0.04 

0.03 
J 0.02 

0.01 




™ I 



I 1 1 iiiii| — I 1 1 n il ] 



I 1 1 iiiii| — I 1 1 n il ] 



cl 



'I' l lli'iliI' M ' l l 



111 

I I I lllll| — I I I lllll| 



III 



-5Slog,(,MS-4.5 

» » « ' 

I I I lllll| — I I I lllll| 



-4Slog,(,MS-3.5 



-I J i i 



I I m ill — I I m ill 



-3Slog,„MS-2.5 



III I 



1.01 0.1 



0.1 



0.1 1 



Fig. 3. — The mass fraction in substructure as a function of 
projected radius {R = R/ R\i) for 3 representative mass bins (rows) 
for low, mean, and high concentration hosts (columns). Each row 
corresponds to the shown mass range of subhalos as a fraction of the 
host halo shown on the far right panel. The y-errorbars correspond 
to the 68 percentile range, and the median of the distribution is 
shown in points. See text for details. 

In Fig. [3] we summarize these effects by showing the 
projected substructure mass fraction as a function of ra- 
dius for the three concentration groups of host halos and 
three representative subhalo mass bins. It is apparent 
that host halos with low concentrations have a higher 
overall projected substructure mass fraction than halos 
with high concentration (see also top panel of Fig. [2]). In 
addition, massive subhalos (subhalos with mass greater 
than 10~* the mass of the host) tend to a median sub- 
structure mass fraction close to zero. However, the upper 
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68 percentile range is increasing relative to the 68 per- 
centile of lower mass subhalos. For example, for masses 
^ 10~^ of the mass of the host, the upper 68 percentile 
implies that up to 16 percent of the low concentration 
halos have substructure mass fractions which are greater 
than 3 percent (with 2% of the low concentration halos 
having a fraction greater than ~ 7.5%). The correspond- 
ing fractions for subhalos with mass ~ 10~^ the mass of 
the host are 1.2 and 3 percent (for the top 16 percentile 
and 2 percentile of hosts respectively) . 

In addition, in high concentration hosts, the 68 per- 
centile range increases with projected radius. The inner 
regions of these halos have evolved for a considerable 
amount of time, while the outer parts reflect more recent 
accretion. On the other hand, in host halos with low 
concentrations, there has been, on average, less time for 
halo substructure to evolve. Moreover, recent accretion 
of subhalos with either highly elongated orbits or with 
orbits that decay rapidly due to dynamical friction can 
deposit large amounts of substructure at either small or 
large halo-centric radii. The consequence of this is that 
the range of substructure mass fractions extends to much 
higher values at all mass bins in low concentration hosts. 

4. STRONG LENSING SIMULATIONS 

We create mock lens systems with characteristics typ- 
ical of observed lens systems, and we include the sub- 
structure populations derived in the semi-analytic cata- 
logs (see previous section). We choose a lensing galaxy 
at z = 0.5 and a source quasar at z = 2 and model the 
lensing halo as a singular isothermal ellipsoid (SIE) with 
projected surface density. 



<0 



-'crit 



2C' 



where ^ is the elliptical coordinate, 



(3) 

^ + y^/l'^i and 

Ecrit is the critical surface density for lensing. In the case 
of a spherical lens, q = 1 and b is the Einstein radius 
of the model. In addition to the smooth lens model, 
an external shear, 7, is applied to account for nearby 
structure, such as a group of galaxies. 9^ is the angle 
between the shear and the major axis of the smooth lens 
model. 

The host halo in the semi-analytic simulations is based 
on a Milky Way-sized dark matter halo. It differs from 
lensing halos in a few respects. First of all, it has no 
baryons. The density profile of the Milky Way in dark 
matter only simulations matches a NFW profile. This 
profile is shallower than isothermal in the inner por- 
tions and steeper than isothermal in the outer portions. 
Real galaxies, however, contain baryons, and observa- 
tions have shown that lensing halos (baryons -I- dark 
matter) follow a roughly isothermal profile near the halo 
center, where stro ng lensing is sensitive to substructure 
( Treu et al.||2006 ). In addition to the differences in den- 
sity profile, we might expect that models with baryons 
would have less substructure than our semi-analytic sim- 
ulations, as the presence of the host galaxy makes tidal 
forces stronger. On the other hand, this effect may be 
small as subhalos that are projected close to the halo 
center may be located at considerably larger physical dis- 
tances. 

The host halo for the substructure catalogs also dif- 



fers with respect to typical lensing halos in a few, more 
subtle ways: the total mass, the formation and merger 
history, and the epoch of observation. The host halo is 
a Milky Way-sized object, while lensing halos are mostly 
those which contain large, elliptical galaxies. We might 
expect that realistic large, lensing halos might have some- 
what more substructure t han Milky Way-sized galaxies 



as larger halos form later |Gao et al. | |2004 
20101). iGao et al.| (|2004|), [Zentner et al.""' 



Giocoli et al. 



WE), and 



van 



den Bosch et al. 1 2005[ ) suggest that the number of suE 



structures scales with host mass as dN/d\nM cx M^ °^, 
which would correspond to only a 15% increase in sub- 
structure compared to our simulations. A similar effect 
is possible as we use substructure catalogs at z = in- 
stead of at z ~ [0.5 — 1] where lensing galaxies are found, 
since larger fract ions of substructur e might be expected 
at earlier epochs (Madau et al. 2008). These effects, how- 
ever, are likely to be significantly small er than th e effect 
of h alo-to-halo variation , as shown by Xu et al. (2009) 
and Zentner et al. (2005). 



Overall, the ettect of baryons on halo structure remains 
an important question, and the optimal method of plac- 
ing dark matter (DM) substructure in a lensing halo re- 
mains unclear. It might be possible to construct a so- 
phisticated model that would connect the properties of 
the host halo to those of the lens halo, but it is nei- 
ther relevant nor necessary for a general study, as is this 
work. Further, treating the host halo and the lensing 
halo rather independently is logical because strong lenses 
have significant surface density contributions from bary- 
onic components of galaxies and because lensing halos 
are described in the data well by isothermal mod els with 
a narrow range of dispersion (Treu et al. 20061). In the 



absence of a model for these effects and all the effects dis- 
cussed above and in order to be conservative in our calcu- 
lations, we preserve (1) the projected distance from the 
center of the halo (i.e., R from the previous section) and 
(2) the fractional mass of substructure from our semi- 
analytic model in our lensing calculations. In addition, 
the host halo mass, Mh = 1.8 x 1O^^M0, is scaled to that 
of a more typical lensing halo, Mions — W^^Mq. Thus, 
the mass of each subhalo is Miens, sub = (Miens /Mh)Msuh 

and aU scale lengths are a;ie„s,sub = (Mions/-Mh)^/^a;sub, 
where x denotes either the tidal radius or the scale ra- 
dius, and both ^lons = R and /ions,sub = /sub- 

Using this model, the Einstein radius of a singular 
isothermal sphere (SIS) lensing halo at z — 0.5 can be 
estimated from Eq. |3j 6 = 0.586". This is generally con- 
sistent with observed lenses which have Einstein radii 



~ 1" (e.g. Browne et al. 2003). It is also a reasonable 
value for any likely lensing halo. The Einstein radius of 
any SIS can be calculated from b — 47r(tT/c)^Dis/Dos, 
where Dig is the angular diameter distance from the 
lens to the source, Z?os is the angular diameter dis- 
tance from the observer to the source, and a is the 1- 
d velocity dispersion. We expect lensing halos to have 
a = [200 - 350] km/s ( |Treu et al.]|2006 ), which corre- 
sponds with b = [0.366 - 1.120J". In Table [l] we show 
the model parameters for the lens halo and environment. 
We refer to these as the macromodel parameters. 

We model subhalos as truncated NFW profiles. There 
is no closed form to calculate the image position deflec- 
tions from NFW profiles with a sharp cut-off at the trun- 
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Name 


b 


<? 


7 


6j (rad) 


Ml 


0.586" 


0.7 


0.2 


7r/3 


M2 


0.586" 


0.7 


0.1 


7r/6 


M3 


0.586" 


0.9 


0.2 


tt/3 



TABLE 1 
Simulation Parameters 



cation radius. We use one of the modified, truncated 



forms provided in Baltz et al. ( 2009 1 , where 



p(r) 



(4) 



where r = rt/rs and rt is the truncation radius. This 
form falls off as fast as beyond the truncation radius 
and the relevant lensing quantities can be calculated sim- 
ply- 

Using the combination of the three different macro- 
models and the set of 200 substructure realizations, we 
perform lensing simulations using an inverse ray-tracing 

code. The position of the source, /3, and the images, 9, 
are related by the lens equation, 



(5) 



where 4> is the lensing potential. In the case of circular 
symmetry in the lens, W<j>u oc M{< R), where M(< R) 
is the mass of the lens enclosed by a cylinder of radius R. 
The lens equation is multi- valued and several image posi- 
tions may correspond to a single source position. Thus, it 
is solved most simply by establishing a grid of image po- 
sitions and calculating the source position for each point 
on the grid. For our simulations, the resolution on the 
image plane is 0.469 mas. 

Every set of three points in the image plane defines 
a triangle and corresponds to a triangle in the source 
plane. The magnification of an image can be calculated 
by the ratio of the area of the image plane triangle to the 
area of the source plane triangle. A common assump- 
tion - one also used in this work - is that the angular 
size of the source is much smaller than that of the dark 
matter substructure. We approximate the source by a 
point. However, for high magnification images, the area 
of the triangle on the source plane may be smaller than 
observed source sizes, as discussed in Sec. [6) 

For each macromodel and each substructure realiza- 
tion, we sample the source plane 50,000 times, creating 
~ 10, 000 mock four image lens systems. We confine our 
investigation to four image systems, since those have a 
sufficient number of observables and have been proved to 
be most useful for substructure studies. 

5. COMPARISONS TO OBSERVATIONAL DATA 

A parameter that can be used to characterize devia- 
tions from a simple, smooth lens is the signed sum of the 
image fluxes from any three adjacent images in a four- 
image system. We work with a scaled version of such a 
sum, 



Rc 



Imi + M2 + Msl 



iMll 



|M3| 



F2 



where fii is the magnification and Fi is the observed flux, 
and both are signed quantities which indicate the image 



1 9 




2 • AO \ 


4 • 


3 y 





Fig. 4. — Sample mock four-image lens, with the lensed images 
(points and labeled by number), the lens center (star), and the 
opening angle for the cusp relation labeled. 

parities. In practice, determining the image parities is 
unnecessary; in a triplet of images, the outer two images 
have one image parity and the middle image has the op- 



posite parity. The ideal cusp relation has Rc 



in 



cases where the lens potential is smooth and the open- 
ing angle A9 spanned by the images from a vertex at 
the lens center approaches zero. An illustrative exam- 
ple of a mock lens system, providing a clear definition 
of A6, is shown in Fig. [4j In practice, we take the sum 
in Eq. ^ over the three images that yield the smallest 
opening angle. The 



20061 |Ai nB,ra et al.||2006| IMetcalf 
Table [5] lists i^cusp values for 



ratio studies (e.g., 


Maccio et al. 


|Metcalf & Amara 


2010). 



lens galaxy and four distinct, point-like images, observed 
in the radio or mid - infrared (noted by "IR" ). This ex- 



cludes £1608-^65 6 (|Koopman s fc FassnachtHl999| |Fass ] 
nacht et al.||2002[ ) BTl27-h3 85 (Koopmfms^'cral. 
and ljl85y-i-154 ( [Myers et al.,.1999^ Kusin et al., 200l| 
which have multiple lensing galaxies! Limiting our- 
selves to distin ct images also excludes MG 2016-1-112 
dGa rrett et al][l 994, 1996; Koopmans et al. 2002), MG 



ehar et al.|1997i|Tonry fc Kochan ck 1999L 
and B1938+666 ( [King et aq[1997| [Tonry fc Kochanek] 
2000| ), in which merging pairs or arcs are seen. We 



specifically include the lens systems which have observ- 
able satellite galaxies near the lens, B2045-I-265 and MG 
0414-1-053, which each have one. Not all of the systems 
in Table [2] are nearly cusp catastrophes - i.e., a cusp 
configuration where A6 — - but in all cases -Rcusp is 
calculated from the 3 images with the smallest opening 
angle. 

Radio data has been traditionally used in lensing stud- 
ies of substructure as it is free of microlensing contami- 
nation from stars in lensing galaxies. This is because the 
size of the quasar source in the radio (~[1 — 10] parsecs) 
is lar ger than the scale size of perturbations from stars 



(6) (e.g., Mao fc Schneider 1998[) and from subsola r mass 



dark matter halos ( [Chen fc Koushiappas 2010). The 



size of the source relative to the perturbations from DM 
substructure, on the other hand, is very small. There- 
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Opening Angle (degrees) 

Fig. 5. — Values of -Rcusp for macromodels Ml (blue dots), M2 
(red dots), and M3 (green dots) without including substructure. 
The observed -Rcusp values from Table [2] are shown as black crosses 
for radio data and as circles for IR data, where mid-IR values are 
used in place of radio values when both are available. 

fore, in our simulations we approximate the source as a 
point. Mid-infrared observations are probably the best 
flux measurements for substructure lensing, as the flux 
measurements are additionally free from extinc tion ef- 
fects, as well as from microlensing effects ( jMinezaki et al 



2004 
have 



|Chiba et al. 2005 1 . Only five lens systems, however, 
jeen observed m the mid-IR (see Table [2]). 



System Name 


A9(°) 


^cusp 


Reference 


1-) 


B2045+265 


34.9 


0.501 


Fassnacht et a 


I.'(!l999i, 










Koopmans et al. (2003) 


2.) 


B0712+472 


76.9 


0.255 


Jackson et al 


(1998), 










Koopmans et al. (2003) 


3.) 


B1422+231 


77.0 


0.187 


Patnaik et al 


(1999), 










Koopmans et al. (2003!) 








0.251 


TChiba et al. { 


2005, IR) 


4.) 


MG 0414+053 


101.5 


0.227 




^^tz et al. 


(1997) 






0.204 


'Minezaki et al. 


(2009 IR) 


5.) 


B1555+375 


102.6 


0.417 


Marlow et al 


(1999), 










Koopmans et 


al. (2003) 


6.) 


B0128-I-437 


123.3 


0.01 


Phillips et al 


. (2000) 


7.) 


PG1115-(-080 


127.5 


0.110 


Chiba et al. ( 


2005, IR) 


8.) 


B1933-I-503 


143.0 


0.39 




Cohn et al. 


(2001) 


9.) 


Q2237-I-030 


146.3 


0.357 




Falco et al. 


T996) 








0.270 




Mir 


lezaki et al.' 


l'f009 IR) 


10.) 


H1413+117 


160.4 


0.22 




MacLeod ct al. 


(2009, IR) 



TABLE 2 
Observed Lens Systems 



Fig. 5] shows the values of i?cusp for macromodels Ml 
(blue dots), M2 (red dots), and M3 (green dots) with- 
out including the effects of substructure. The observed 
i?cusp values from Table [21 are shown as black crosses and 
circles, where mid-IR values are used in place of radio 
values where both are available. Small opening angles 
correspond to small cusp relation values, as is expected. 
At larger opening angles, for each model the values of 
-Rcusp bifurcate, corresponding to source positions along 
the minor axis (smaller i?cusp values) or the major axis 
(larger i?cusp values) of the lensing halo. As expected, 
values of i?cusp dep end on the macromodel parameters 
( [Keeton et al.| 2003). At the smallest of opening angles, 
-n^cusp 0, regardless of macromodel. At opening an- 



gles, < 90°, the observed values of -Rcusp clearly deviate 
from what may be expected from smooth models. 

The effect of substructure on the cusp relation is shown 
in Fig. [6j The no substructure results from Fig. [5] are 
reproduced here in red and orange shaded areas. The 
borders of the colored areas are contours of equal prob- 
ability, and the red and orange patches encompass the 
areas where 99% and 100% of the cusp relation values 
are found. These areas are well contained within a small 
portion of the -Rcusp — parameter space. Including 
the effect of substructure gives values which are shown 
by the dark gray and light gray shaded areas. The dark 
gray area, which encompasses the area where 99% of the 
values are found, overlaps significantly with the region 
in which the no substructure values lie. And where the 
substructure-included values do extend beyond this area, 
they mostly lie fairly close to the no substructure values. 
In these cases, the deviations from the no substructure 
case are not due to small perturbations to individual im- 
ages but perturbations that effectively change the macro- 
model of the lens systems - e.g., shifting the lens center 
or changing the overall ellipticity of the lens. The light 
gray region shows where the final 1% of -Rcusp values lie 
for the substructure-included systems. Here, some large 
deviations from the no substructure case can be seen. In 
particular, some small number of lenses show -Rcusp ^ 1 
at opening angles less than 90° . So large deviations from 
the cusp relation due to the influence of substructure are 
possible but are infrequent. 

We bin the cusp relation values by opening angle and 
focus on three bins that contain observed cusp systems; 
30° < A6I < 40°, 70° < A9 < 80°, and 100° < AO < 
110°. When comparing our mock lens systems to ob- 
served -Rcusp values, we find that the probability of high 
-Rcusp values remains small. A summary is presented in 
Table|3l Here we use all 200 substructure realizations and 
estimate the cumulative probability of observing -Rcusp 
values that are greater than provided by smooth models, 
-P(^cusp > RS), where RS is the -Rcusp value in the no 
substructure case where P(-Rcusp > -R5'|no subhalos) = 
0. For mock lens systems using the substructure realiza- 
tions, P(-Rcusp > RS) ~ 10^^. In addition, the cumula- 
tive probability of observing -Rcusp values that are greater 
than the observed values, P(-Rcusp > RO) is estimated, 
where RO is the largest -Rcusp value observed for a given 
bin of opening angle. For mock lens systems using the 
substructure realizations, P(-Rcusp > -RO) ~ 10^^. 

In Fig. [Tj we show the distribution of -Rcusp values in 
a given bm of opening angle. Here, the distribution of 
values for mock lenses with no substructure cluster near 
small -Rcusp values (shown as cyan bars). On the other 
hand, the observed values are found at large -Rcusp values 
(cross-hatched, magenta bars). The distribution of val- 
ues for mock lenses with substructure are shown in red 
bars and black-lined bars. Substructure induces high- 
-Rcusp tails which extend to (and, in some cases, past) 
the observed values of the cusp relation. In all three 
opening angle bins, these tails contain only a small frac- 
tion of the lens systems (a few percent). Interestingly, 
the high- -Rcusp tails extend to larger values of -Rcusp for 
smaller opening angles than for larger opening angles. 
The substructure realizations are split by the concentra- 
tion of the host galaxy, as discussed in Sec. [2l plotting 
only the realizations that constitute the srnallest 16% 
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A9 [degrees] AS [degrees] AS [degrees] 

Fig. 6. — Values of i?cusp for Models 1 (left panel), 2 (center panel), and 3 (right panel) including all substructure realizations (light and 
dark gray areas) and overplotted with no substructure realizations (orange and red areas). Red and dark gray areas show where 99% of 
the values lie, while light gray and orange show where the remaining 1% of the values are found. The observed -Rcusp values from Table[2] 
are shown in blue dots for radio data and green asterisks for mid-IR, where mid-IR values are used in place of radio values when both are 
available. 





Macromodcl 


^cusp (R'S 


a 




P(i?cusp > RS) 


Rcusp (ROf 




P(Rcusp > RO) 


30 - 40° 


Ml 


0.1 


0.110 


0.5 


0.042 




M2 




0.048 




0.001 




M3 




0.042 




0.001 


70 - 80° 


Ml 


0.15 


0.0198 


0.25 


0.0121 




M2 




0.0320 




0.0138 




M3 




0.0148 




0.0062 


100 - 110° 


Ml 


0.25 


0.0070 


0.4 


0.0008 




M2 




0.0171 




0.0028 




M3 




0.0052 




0.0014 



^The largest Rcusp value in models with no substructure. 
''The largest observed value of Rcusp- 

TABLE 3 

Cumulative Probability of Cusp Relation Values 



of host concentrations (red bars) and the reahzations 
that make up the largest 16% (black-hned bars). We 
see that the correlation between host concentration and 
amount of substructure persists into observable lensing 
results. The hosts with the smallest concentrations have 
the largest amount of substructure and larger values of 
i?cusp when compared to the hosts with the largest con- 
centrations. This is best seen in the bins of 70 — 80° and 
100—110°, where the number of mock lens systems num- 
ber in the thousands. In the smallest opening angle bin 
of 30 — 40° , the number of mock lens systems is less than 
20 per substructure realization and differences based on 
halo concentration are more difficult to discern. 

We test which parameterizations of the substructure 
content of lensing galaxies are most sensitive to anoma- 
lous fluxes and are most useful for comparing observed 
lens systems to the results of simulations. First, we ex- 
plicitly tailor the definition of /sub from Section [2j Eq. [T] 
to match our macromodcl conditions: 



(7) 



where Mi{Ri < 2b) is a subhalo with mass M that is 
found to be projected within twice the Einstein radius of 
the lens b (within which all the images using macromod- 
els Ml, M2, and M3 have been found) and Mio„s(i? < 26) 
is the mass of the lensing galaxy within 26. In this equa- 
tion, we have the freedom to define what it means for 
a subhalo to be projected within R — 26. As discussed 



previously, the substructure in our catalogs have density 
profiles which are described by a truncated NEW pro- 
file, which has three parameters: the concentration, the 
bound mass, and the tidal radius. We could chose to 
include any subhalo which is projected such that some 
portion of the area within R = 2b falls within the tidal ra- 
di us of the subhalo. This parameterization is supported 
by Rozo et al. ( 2006 ) , who suggest that the most impor- 
tant UM halos for niagnification perturbations are those 
within a tidal radius of an image. Another possibility is 
that a single subhalo could be projected such that some 
portion of the area within R — 2b falls within the scale 
radius of the subhalo. A final possibility is that a single 
subhalo could be projected such that the center of the 
subhalo falls within the area defined by i? = 26. This is 
the parameterization used in previous studies, when ana- 
lyzin g the outpu ts of dark matter simulations, ( Mao et al 
2004;iXu et al.f2009). The fist of different parameteriza- 
tions, then, is^as follow: 1.) include subhalos whose tidal 
radii fall within 26; 2.) include subhalos whose scale radii 
fall within that area; and 3.) include subhalos whose halo 
centers fall within that area. 

In Eig. Isj the substructure realizations are split into 
low /sub and high /sub bins using the three criteria. Note 
that some realizations have no subhalos which satisfy the 
criterion, i.e., /sub = 0: 4 out of 200 for case 1, 45 out of 
200 for case 2, and 58 out of 200 for case 3. In addition, 
note that values of /sub > 1 are possible. This may be 
the result of two effects: MicnsiR < 26) from Eq. [Tjonly 
refers to the smooth mass of the host halo and not the 
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tions where the host halo concentration is in the lowest 16 percent 
(red bars) and in the highest 16 percent (black lines). At each 
Kcusp bin, the bars show the results for Ml, M2, and M3 from 
left to right. The cyan colored bars show the results using the M2 
macromodel parameters and no substructure. The cross-hatched 
magenta shows the observed results. From top to bottom, the 
panels show opening angles of 30-40°, 70-80°, and 100-110°. 

total mass, and we use ttie full subhalo mass even when 
most of the subhalo's bound mass falls outside of the rel- 
evant region. Using all three criteria, we can see that 
substructure realizations with high /sub values also have 
larger and more frequent violations of the cusp relation, 
when compared to low /sub realizations. The differences 



between low and high values of /sub are greatest for val- 
ues of /sub calculated via halo centers (case 3). Here, 
the median /sub value is /sub = 0.13%, and substructure 
realizations with /sub greater than the median have an 
appreciable but small 10~^) fraction of large i?cusp 
values while those with /sub smaller than the median 
may have no large cusp relation values. 

The value of the cusp relation for every mock lens sys- 
tem is dependent on the relative positions of the subhalos 
and the source. The placement of the subhalos and the 
source is stochastic ~ depending upon the merger history 
of the lens halo, the direction in which the halo is pro- 
jected, and the sampling of the source plane. Thus, the 
results of Fig. [8] could be due to just a few substructure 
realizations that always give high-i?cusp values, while the 
majority of substructure realization give small i?cusp val- 
ues regardless of /sub • We address this possibility in Fig. 
[9j For every bin of /subi we measure the fraction of sub- 
structure realizations that have at least one large cusp 
relation value among their set of mock lenses. A small 
fraction of the realizations with small values of /sub have 
large i?cusp values, while a large fraction of the realiza- 
tions with large values of /sub have large i?cusp values. 
This behavior is most pronounced for a /sub definition 
that includes only subhalos whose center falls near the 
image positions, even when considering that, using this 
definition, the bins with /sub > 0.05 contain only a few 
realizations each. Thus, /sub is reliably correlated with 
i?cusp values. 

The parameterization of substructure in lensing halos 
which is most often used in studies of anomalous flux 
ratios is /sub- Other possibilities for parameterizing sub- 
structure exist: the number of subhalos and the aver age 
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subhalo massj^ We test these possibilities in Fig 
Here, subhalos whose halo centers fall within R — 2b are 
counted and the mean subhalo mass measured. The sub- 
structure realizations with a larger number of subhalos 
also have larger and more frequent violations of the cusp 
relation, when compared to realizations with a smaller 
number of subhalos. In addition, the substructure real- 
izations with a larger average subhalo mass have larger 
and more frequent violations of the cusp relation, when 
compared to realizations with a smaller average subhalo 
mass. However, when compared to Fig. [Sj neither the 
number of subhalos or the mean subhalo mass is as dis- 
criminating as /sub- For example, for the observed value 
of i?cusp — 0.25 (for opening angle bin 70 — 80°), Figure 
|10| suggests that the mean subhalo mass and the num- 
ber of subhalos could be large or small. On the other 
hand, Figure[8]is clear in showing that Luh must be large. 
This r esult is consistent with that of IMetcalf fc Amaral 
(20101, who find that the frequency of flux anomalies is 
a function of /sub and not of the number or average mass 
separately. 

6. CAVEATS 

In this section we discuss some specific caveats to our 
calculation that we are aware of and should be taken into 
account in the interpretation of these results. 

* Note that the average subhalo mass is not uniquely defined be- 
cause it scales with the limits of the subhalo mass function, i.e., the 
average mass can be changed by changing the lower limit of subhalo 
masses investigated. The lower limit in our work is 10~^Mip„„. 
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Fig. 8. — Top: Distribution of f^uh values for three different criteria. See text for additional details. The black vertical lines show the 
sample median. Some substructure realizations have f^uh = and are not shown in the panels. Bottom: The corresponding distribution of 
ficusp values for AO = 70 — 80°. The substructure realizations where f^ub is greater than the median are shown with filled red bars, and 
realizations where /sub is less than the median are shown with black lines. In each Rcusp bin, the results from Ml, M2, and M3 are shown 
from left to right. 



In Sec.|4| we discussed how the simulations of substruc- 
ture differ between the semi-analytic host halo and the 
typical lensing halo and the accommodations we make in 
order to incorporate the substructure catalogs into our 
mock lenses. We make simple adjustments ~ adopting the 
projected distance distribution and projected mass frac- 
tion from the simulations while scaling the mass of the 
host galaxy from Mh = 1.8x W'^Mq to Mio„s = 10"Mq. 
The scaleability of the halo and subhalo masses is in line 
with the results of dark matter num erical simulations 
(iGao et al.||2004| iGiocoh et al.||2010|). The host halo 



in the semi-analytic realizations also differs with respect 
to typical lensing halos in more subtle ways: the for- 
mation and merger h istory and the epoch of observation 
(Madau et al. 2008). While these differences are not 
likely to significantly effect the dist ance distributio n or 
mass fraction of substructure (e.g., Xu et al. 2009 and 
as discussed in Sec. they have unquantified effects on 
the internal structure of the subhalos - their scale and 
tidal radii. The concentration of a subhalo, however, is a 
slowly varying function of mass and time, so that biases 
in this analysis induced by these effects are likely to be 
small compared to natural halo-to-halo variation. 

With regards to finite source effects, in this work we use 
the point source approximation, which - as stated pre- 
viously - could bias the calculated m agnification of high 
magnification images. For example, [Metcalf fc Amara] 
(12010) suggest that finite source effects could have large 
efi'ects on the magnifications of images. For lens sys- 
tems with substructure and total magnifications above 
fi ~ 10, the fractional deviation in magnification for a 
point source and for a 1 pc source approaches 1. For 
smooth lens halos, the same holds true for magnifica- 
tions above fi = 20. The bias in magnifications most 
often manifests itself as an overestimate of the magnifi- 
cation of an image when a point-like source is used in 



place of a finite source. 

However, tests using our code fail to find significant 
biases due to the point-like source assumption. This dif- 
ference can be accounted for by two effects. First of all, 
we consider only the cusp relation, which is the difference 
between magnifications of images, mi tigating the effects 
of any finite source effect. Secondly, [Metcalf fc Amar"a 
(20101 use a model of substructure where the surface 
density profile is a power-law E cx and a = 1 or 

a = 0.5. This profile is significantly steeper than the 
NFW profile used in this work and is likely to cause big- 
ger magnification perturbations than the subhalos in our 
work. 

Finally, with regards to magnification bias, the dis- 
tribution of opening angle values in our mock lenses, as 
shown in Figs. [5]and[6] is not representative of the distri- 
bution expected in observed lens samples. For example, 
our mock data set has fewer lenses with small opening 
angles than expected observationally. This is because the 
distribution of opening angles in our mock samples is a 
result of a combination of the macromodel parameters 
and the manner in which the source plane is sampled in 
order to create the mock lens systems. We sample the 
source plane uniformly, a method which does not account 
for the magnification bias of observations - the fact that 
bri ghter systems are more like ly to be observed. Follow- 



ing peetorT^^^ZaHu^ (2004), it can be shown that the 
probability density in tne source plane is related to the 
luminosity function of sources, such that 



p{u)du = n{uy ^du. 



(8) 



when the luminosity function of sources has a power-law 
form dN/dS cx 5'"'^, where S is the source flux and where 
u is the position of the source. Thus if = 2, the prob- 
ability of observing a lens system with source position u 
is proportional to the total magnification of the system. 
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Fig. 9. — Fraction of realizations in each bin of /sub which also 
contain -Rcuap > ^?cut values. Three different criteria for /sub are 
considered: including substructures whose tidal radii fall within 
R = 2b (black bars); including substructures whose scale radius 
fall within R = 2b (red bars); and including substructures whose 
halo centers fall within R = 2b (blue bars). Opening angles of 
70-80° are shown and results using macromodels Ml, M2, and 
M3 are combined. The threshold values, Rcut = 0.15 (top) and 
Rent = 0.25 (bottom) are taken from Table p] and represent the 
ficusp value that is greater than provided by smooth models and 
ficusp value that is greater than the observed values, respectively. 



At. Stu dies of the Cosmic Lens All-Sky Survey (C LASS; 
Browne et al. 2003 ) have sug gested that i/ = 2.1 (|Rusin| 
fe I'egmark 2001; C hae|2003| . The total magnificafionof 
a system is related to the opening angle of a cusp so that 
the smallest opening angles are found in the systems with 
the largest magnifications. In our observational sample, 
one of ten in the observational set has an opening angle 
less than 40°. The mock lens systems, in contrast, have 
very few systems (< 1%) with opening angles less than 
40°. Accounting for magnification bias would increase 
the number of systems with opening angles less than 40° 
to a few percent. Given the effect of magnification bias, 
we negate its effect by making comparisons of observed 
lens systems to our mock lens systems using only data 
that is binned in opening angles. 

7. CONCLUSIONS AND DISCUSSION 

Dark matter substructure is a generic prediction of the 
cold dark matter model and probes of substructure are a 



key test of the nature of the dark matter. Strong gravi- 
tational lensing can be used to detect and measure dark 
matter subhalos, as they cause perturbations to the im- 
age magnifications, positions, and time delays in lens sys- 
tems with galaxy lenses. Putting robust constraints on 
these observational results require robust predictions for 
CDM substructure that account for halo-to-halo vari- 
ation. Current predictions from numerical simulations 
have used a limited number of simulated halos. We in- 
vestigate the question of halo-to-halo variation in sub- 
structure using semi-analytical models and test how it 
affects perturbations to the magnifications of images, as 
measured by the cusp relation parameter, i?cusp- We find 
the following results. 

1. The total projected substructure mass fraction for 
halos with mass greater than 10~^ of the mass of 
the host is ^ 9% with a 95 percentile in the range 
[2 — 17]%. Low concentration halos have a higher 
projected mass fraction 20%), with a signifi- 
cantly larger 95 percentile range. High concentra- 
tion halos have a lower mass fraction (~ 5%), and 
a much smaller 95 percentile range. 

2. We find that the mass fraction within the inner 
3% of the projected radius is fsnh = 0.25%, with 
a 95 percentile in the range of fsnh = [0—1]% for 
the complete sample of 200 realizations. For low- 
concentration host halos the median substructure 
mass fraction is /sub ^ 0.5% with a 68 percentile 
in the range of fsuh — 0.2 — 0.8%. These results are 
genera lly consistent with the results of |Xu et al.| 
| 2009| , who find .01% < f^uh < -7%, using 6 ha- 
los from the Aquarius simulation, as well as with 
all previously published limits on /sub- The mass 
fraction of host halos with higher concentrations is 
lower, with /sub < 0.1% for host halos with con- 
centrations in the upper 16% of the distribution. 

3. Accounting for halo-to-halo variation, the observed 
values of i?cusp are greater than the values pre- 
dicted by smooth lens models, and larger than 
99% of the values predicted by substructure mod- 
els. Given halo-to-halo variation, the probability 
of observing cusp relation values greater than sup- 
ported by smooth models is ~10~^ and the proba- 
bility of observing cusp relation values greater than 
the largest currently observed values is ^10""^. 

4. Substructure fractions are correlated with the con- 
centrations of the host halos. Lower-concentration 
host halos have larger substructure fractions so 
that lower concentration host halos give rise to 
more substructure and more violations of the cusp 
relation. 

5. For lensing studies of substructure, the substruc- 
ture content of lensing halos is best parameterized 
by a definition of the substructure mass fraction, 
/sub, which includes only the subhalos whose halo 
centers fall within the area in which images are ex- 
pected to be found (twice the Einstein radius in 
our models, ^ 3% of the virial radius of the lens 
halo, generically) and which excludes more distant 
subhalos that fall within a scale or tidal radius of 
a lensed image. 



12 



Chen, Koushiappas & Zentner 



0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




1 23456789 
Number of Substructures 



10 




1e-05 0.0001 0.001 0.01 
Mean Substructure Mass (/M^oJ 




Fig. 10. — Top: The number of subhalos whose centers fall within R = 2b (left), and the mean mass for those subhalos (right). The 
black vertical lines show the sample median. Bottom: The corresponding distribution of iJcusp values for A9 = 70 — 80° . The substructure 
realizations whore number of mass is greater than the median is shown in red bars, and realizations where it is less than the median in 
shown in black lines. In each -Rcusp bin, the results from Ml, M2, and M3 arc shown from left to right. 



6. /sub is a better predictor of -Rcusp when compared 
to the number of subhalos or the mean subhalo 
mass. 

We attempt to isolate the effect of halo-to-halo vari- 
ation in substructure populations from macromodel ef- 
fects and, thus, do not model any particular observed 
lens system. Our efforts can still be compared qualita- 
tively to previous attempts to constrain the substructure 
mass fraction. For example, Dalai & Kochanek (20021 
find that /sub = 0.02 with b% < /sub < 7% -^Mmi 
90% confidence intervals. This result is marginally con- 
sistent with predict ions of /sub in previous works and 
in this paper, while Metcalf & Amara (2010) find rough 
consistency with such predictions. Both of these works 
use substructure models with density profiles that are 
steeper than the NFW profile used in this work. Steeper 
profiles are more effective at causing magnification per- 
turbations. Thus, it is possible that both these works 
are underestimating the /sub that is supported by lens- 
ing data. Other works have also suggested that lensing 



(IMaccio et al. 


2006 


Amara et al. 


Xu et al. ( 


2om 


suggest that, using 



2006 |Chen||2009[ ) _ 

the 6 haios in the Aquarius simulation, the probability 



In summary, we find that when accounting for halo-to- 
halo variation the probability of observing i?cusp values 
that match observations may be as large as ~10~^, but 
the tension between observations and numerical simula- 
tions persist. This tension may be eased by assuming 
larger values of the substructure mass fraction. Possible 
solutions, then, are 1.) additional substructures which 
are unresolved in current numerical and semi-analytic 
simulations and 2.) biases which cause the observational 
sample to prefer lensing halos with high substructure 
mass fraction or, relatedly, low concentration. Both these 
possibilities, however, are speculative. 

While this work only probes the effect of subhalos with 
masses ~1O^M0 or larger, others h ave studied smaller 
substr uctures. Using a ACDM model, 
(20061 investigated the role of sub' 
lO^lO^Me 

ing large cus p relation violations. The Aquarius simula- 
tion used by Xu et al. ( 2009) resolves subhalos down to 



Maccio & Miranda 
alos 



lalos with masses 
and found them to be insufficient in creat- 



of the observed values of Rn 



^10- 



10°h~ Mq. The cumulative effect of smaller and larger 
subhalos produces their result of a probability that ob- 
servations match expectations from ACDM of ~ 10~^. 

Currently, there are no known observational biases 
which would prefer high substructure or low concentra- 
tion lensing halos. If the observational bias from the con- 
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centration of lensing lialos is conside red alone, Iiigh con 



centr ation halos should be preferred ( Mandelbaum et al. 
20091. In addition, the relation between observational 
biases, substructure content, and other characteristics of 
lenses may be complicated. For example, triaxial halos 
projected along the long axis are better strong lenses in 
general and may have more substructure projected near 
images ( Zentner||2006 1 , but triaxial halos along the mid- 



dle axis ar e better at creating four- image lenses (Rozo 
et al.||200"7 ). Accounting for all the biases due to lensmg 
halo characteristics is difficult and requires future inves- 
tigation. 
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Brown University. ARZ is funded by the National Sci- 
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